Diffusion entropy analysis on the stride interval fluctuation of human gait 
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In this paper, the diffusion entropy technique is applied to investigate the scaling behavior of 
stride interval fluctuations of human gait. The scaling behavior of the stride interval of human 
walking at normal, slow and fast rate are similar; with the scale-invariance exponents in the interval 
[0.663, 0.955], of which the mean value is 0.821 ± 0.011. Dynamical analysis of these stride interval 
fluctuations reveals a self-similar pattern: Fluctuation at one time scale are statistically similar to 
those at multiple other time scales, at least over hundreds of steps, while the healthy subjects walk 
at their normal rate. The long-range correlations are observed during the spontaneous walking after 
the removal of the trend in the time series with Fourier filter. These findings uncover that the fractal 
dynamics of stride interval of human gait are normally intrinsic to the locomotor systems. 

PACS numbers: 87.90.+y,89.75.Da,05.40-a,05.45.Tp 



I. INTRODUCTION 

Recently, it has been recognized that in many natu- 
ral sequences the elements are not positioned randomly, 
but exhibit long-range correlations and fractal dynam- 
ics. Prominent examples include noncoding DNA se- 
quences , human heartbeat 0, 01 j human brain elec- 
troencephalogram pL teenagers' births in Texas p|, and 
financial time series [6j . The common feature of all these 
diverse systems with long-range correlations is that scal- 
ing behavior decay by a power law, where a characteristic 
scale is absent, and the scaling behavior is usually rather 
robust and universal. 

The scaling behavior in complex systems is not only in- 
teresting in physical sense, but also provides an intrinsic 
description of the system, and highlights the system dy- 
namical mechanism. Several useful variance-based meth- 
ods, such as the probability moment method and the 
fluctuation approach as well as the de-trended fluctuation 
approach pj , are proposed to detect the scale-invariance 
properties. However, these variance-based methods have 
two basic shortcomings. One is that the scale-invariance 
property can be detected but the value of the exponent 
cannot be obtained accurately. The other is that for some 
processes, like the Levy flight, the variance tends to infi- 
nite and these methods are unavailable at all H, @ • Al- 
though the infinite can not be reached due to the finite 
records of empirical data, clearly we will obtain mislead- 
ing information about the dynamics under these condi- 
tions. 

The dynamics of human gait is relative to locomotor 
system's synthesizing inputs from the motor cortex, the 
basal ganglia and the cerebellum, as well as feedback 
from the vestibular, visual and proprioceptive sources. 



Therefore, to investigate the scaling behavior of the stride 
interval fluctuations of human gait is very interesting 
with both the physical and physiological communities. 
By using some recently proposed approaches for non- 
linear data, some scientists studied the statistical char- 
acters of the stride interval fluctuations of human gait 
E3, EH O O [3 EE HE HI Hausdorff et al [H ]n| 
demonstrated that the stride interval time series exhibits 
long-range correlations, suggesting that the human gait 
displays a self-similar activity. Subsequent studies, by 
West et al and Goldberger et al [H El Q 13, sup- 
ported the above conclusion, and several dynamics mod- 
els arc established to mimic the human gait pH , fli^ . 
Furthermore, Perc [nj found that the human gait pro- 
cess possesses some typical properties like a deterministic 
chaotic system. In this paper, we apply the diffusion en- 
tropy (DE) technique to accurately detect and obtain the 
scaling property of stride interval fluctuations of human 
gait on a Fourier analysis. 



II. DIFFUSION ENTROPY TECHNIQUE 
BASED ON FOURIER ANALYSIS 



To overcome the mentioned shortcomings in the 
variance-based methods, Grigolini et al. @, Q designed 
the approach named diffusion entropy analysis (DEA). 
To keep our description as self-contained as possible, we 
briefly review the DEA method. 

Consider a complex system containing a large amount 
of particles. The scale-invariance property in the diffu- 
sion process of this system can be described mathemati- 
cally with the probability distribution function as 



P(x,t) 



(1) 
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where x is the displacements of the particles in the com- 
plex system and S the scale-invariance exponent. The 
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FIG. 1: A typical example about the stride interval fluctua- 
tion of human gait. 



theoretical foundation of this property is the Central 
Limit Theorem and the Generalized Central Limit Theo- 
rem [2(], EI ■ For 5 — 0.5, the diffusion process is the stan- 
dard diffusion and F(y) is the Gaussian function. And 
5 7^ 0.5 indicates the deviation of the dynamical process 
from the normal stochastic one. For example, the case 
with 5 > 0.5 corresponds to the diffusion process of Levy 
walk. 

For a time series, the delay-register vectors, denoted 
with {y k ,y k+ i, ■ ■ ■,y fe+m _ 1 |A; = 1, 2, 3, • • •, N - m + 1}, 
can be regarded as the trajectories of N — m + 1 particles 
during the period of to m. fn this way we can map 
a time series to a diffusion process, called overlapping 
diffusion process here. An alternative solution is to sepa- 
rate the considered time series into many non-overlapping 
segments and regard these segments as the trajectories. 

To make the data suitable for the illustration of scaling 
behavior at all scales, we propose a method based on 
the Fourier transform to detrend the nonstationary trend 
in whole temporal domain. First, transform the time 
series to Fourier space and then, use the invert Fourier 
transform to obtain the time series with cutoff coefficient. 
We consider the curtailed time series as the trend and 
subtract it from original series. Then, the final time series 
is regarded as a steady series, whose overlapping diffusion 
process reads 

k+t 

a*(t)=X;< J -,*=l,2,.-.,W-t+l. (2) 

j=k 



Consequently, the Shannon entropy can be defined as 

/+oo 
P(x,t)log 10 [P(x,t)]dx. (3) 
-DO 

A simple algebraic leads to 

S(t)=A + Slog w (t), (4) 

where 

/+oo 
F(y)]og w [F(y)]dx 1 y=- r (5) 

The DEA method has been used to analysis many time 
series in different research fields, such as the solar in- 
duced atmosphere temperature series l22l| . the intermit- 
tcncy time scries in fluid turbulence 23], the spectra of 
complex networks |2 11 , the output spike trains of neurons 
|2^j . the index of financial market [2£j, and so on. 

III. DATA ANALYSIS 

Herein, we map a time series to a diffusive process and 
introduce the DE technique to investigate the time series 
of stride interval fluctuations of human gait, which is ob- 
tained from healthy subjects who walked for I hour at 
their normal, slow and fast paces |27j . The data contains 
the stride interval fluctuations of ten young healthy men, 
given an arbitrary ID (siOl, sz02, si03, • ■ • , silO). Partici- 
pants have no history of any neuromuscular, respiratory, 
or cardiovascular disorders and are taking no medica- 
tion. Mean age is 21.7 yr (range 18 — 29 yr). Height is 
1.77±0.08(SD)m, and weight is 71.8 ± 10.7(SD)fc#. All 
the subjects provided informed written consent. Subjects 
walked continuously on level ground around an obstacle 
free, long (either 225 or 400 meters), approximately oval 
path and the stride interval is measured using ultra-thin, 
force sensitive switches taped inside one shoe. A typical 
example is shown in Fig. 1. 

The stride interval significantly decreased and the ve- 
locity significantly increased with each chosen walking 
rate, as expected. The mean stride intervals of the three 
subjects are 1.3 ± 0.2s, 1.1 ± 0.1s, and 1.0 ± 0.1s during 
the slow, normal, fast walking trials, respectively. And 
the mean velocities are 1.0 ± 0.2to/s, 1.4 ± O.lm/s, and 
1.7 ± 0.1m/ s during the slow, normal, and fast walking, 
respectively. The mean velocity increases by an average 
of 77% from slow to fast walking. A wide range of walk- 
ing rates are obtained, enable us to test for the effects 
of walking rate on the scaling behavior indicating the 
long-range correlations. 

The locomotor control system maintains the stride in- 
terval at an almost constant level throughout the 1 hour 
walking. Nevertheless, the stride interval fluctuations is 
in a highly complex, seemingly random fashion. In order 
to truly uncover the scaling behavior of the stride interval 
fluctuation of human gait, we study all the 30 samples 
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FIG. 2: (Color online) The scaling behavior of the stride 
interval fluctuations of human gait. (A, (B) and (C) denote 
the the subject Fast, Normal and Slow respectively. 



classifying into three subjects: Fast, normal, and slow. 
The scaling behaviors of the three classes of subjects, 
obtained by using DE analysis, are shown in Fig. 2A, 2B 
and 2C, respectively. Each figure contains 10 samples, 
of which the time scale is from 10 to 300. The results 
indicate that the stride interval time series is not com- 
pletely random (uncorrelated) , instead, it exhibits the 
scale-invariance property and long-range correlation at 
all the three walking rates. 

Furthermore, Fig. 3 illustrates the dependence of S on 
self-determined walking rate for all three subjects. For 



c 
a) 
c 
o 

I 

CO 

c 

w 
u 

V) 

2 

Q 




Upward:Original time series 
Downward:Randomly Shuffled 




Normal Slow 

Walking condition 



FIG. 3: (Color online) Dependence of S on self-selected walk- 
ing rate. 



the 29 samples of 30 1-h trials, except one samples si08 
of the slow subject (see Fig.2C), the scale- invariant ex- 
ponents S is around 0.821 ± 0.011 (range 0.663 to 0.955). 
Thus for all subjects at all rates, the stride interval time 
series displayed scaling behavior and long-range power 
law correlations. 

There are a strong similarity of S on chosen walking 
rate: S is 0.829±0.011, 0.827±0.011, 0.809±0.013 during 
the fast-, normal-, slow-walking trails, respectively (see 
Fig. 3). In principle, the estimation of 6 depends on the 
fitting range (il,i2). It is preferable to have t% 3> t\. 
However, due to the finite size of data, a larger ti may 
lead to a bigger error in 8 28]. In this paper, all the 
results are obtained by linear fitting within the region 
(10,300). 



IV. CONCLUSION 

In summary, by means of the DE method we inves- 
tigate the scaling behavior embedded in the time se- 
ries of stride interval fluctuation. Scale-invariance ex- 
ponents of the three subjects are almost the same, being 
in the interval of [0.663, 0.955], of which the mean value 
is 0.821 ±0.011. Dynamical analysis of these step-to-step 
fluctuations reveals a self-similar pattern: Fluctuation at 
one time scale are statistically similar to those at mul- 
tiple other time scales, at least over hundreds of steps, 
while the healthy subjects walk at their normal rate. The 
long-rang correlation is observed during the spontaneous 
walking by removal of the drift or trend in the time series. 
Thus the above features uncover the fractal dynamics of 
spontaneous stride interval are normally intrinsic to the 
locomotor system. 
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